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Abstract 

To verify the conjecture that Yang- Mills theory in the infrared limit is equivalent to 
a spin system whose excitations are knot solitons, a numerical algorithm based on the 
inverse Monte Carlo method is proposed. To investigate the stability of the effective 
spin field action, numerical studies of the renormalization group flow for the coupling 
constants are suggested. A universality of the effective spin field action is also discussed. 

1. It was conjectured [lj that the SU(2) Yang-Mills theory in the infrared region can be 
described as a spin system with the following action 

S = Jdx {m 2 (^n) 2 + g- 2 [n ■ (c^n x d u n)] 2 } , (1) 

where n 2 = 1 (boldface letters stand for three- vectors) . The path integral representation 
of the effective action (|I]) can be deduced from the Yang-Mills theory path integral via an 
implicit change of integration variables @. The analysis can be extended to the SU(N) 
case [^|, |J. Nonperturbative excitations of the effective theory are knot solitons Knot 
solitons look more like stringy excitations, which is believed to be a right physical picture 
of nonperturbative excitations of gauge fields. Yet, if the effective action ([!]) turns out to 
be a good approximation to the Yang- Mills theory in the infrared limit, the nonperturbative 
dynamics can be studied by quantum soliton theory methods. The mass gap in the spectrum 
of quantum Yang-Mills theory would therefore naturally be introduced as the lowest energy 
bound in the quantum soliton spectrum. A few important questions are to be addressed to 
validate or invalidate this attractive picture. First, what is the actual value of the mass scale 
m 2 which determines the low energy bound (in the classical soliton theory) ||? Second, is the 
effective action stable from the point of view of the renormalization group flow of its coupling 
constants? Third, how big are the higher order corrections to (|l|)? The purpose of this paper 
is to set up a numerical approach to answer these questions. 

The existence of degrees of freedom whose dynamics dominates in the infrared region 
of Yang-Mills theory was established in numerical simulations 0, [8|, || of lattice Yang-Mills 
theories some time ago. It was observed that dominant contributions to the string tension 
come from topological defects (monopoles) which occur in typical vacuum configurations of 
gauge fields when the latter are taken in a special gauge known as the maximal Abelian gauge. 
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Topological defects unavoidably occur in any gauge that breaks the gauge group to its maximal 
Abelian subgroup [|K|. In classical theory, it is evident from the fact that the homotopies of 
the gauge group and its maximal Abelian subgroup are different (see more on the gauge 
fixing problem in quantum gauge theories in [0). The importance of the above numerical 
discovery is that the defects alone are sufficient to reproduce essential nonperturbative features 
of Yang-Mills theory. The numerical procedure of singling out topological defects in Abelian 
projections of the lattice Yang-Mills theory can therefore be used in a theoretical analysis to 
parameterize the relevant degrees of freedom of the Yang-Mills connection and compute their 
effective action. It should be noted that the gauge fixing here does not serve its conventional 
purpose - removing nonphysical degrees of freedom - but rather it becomes an auxiliary tool 
to identify the degrees of freedom relevant for the infrared physics of the Yang-Mills theory. 
After a change of variables, that splits all the Yang-Mills degrees of freedom into the "infrared 
relevant" ones and the rest, is found, the effective action can be computed in any convenient 
gauge and its gauge invariance can be established by the standard BRST technique || [|. 

In our earlier works [Q, [J, a relation between the spin field n and topological defects of 
the connection in the maximal Abelian gauge has been found 

A M = g^dpn x n + nCf, + , (2) 

where W M satisfies the following conditions. It is perpendicular to n and 

d^W^ + + nC;) x = V M (« + nC)W ^ = . (3) 



Here ct M = g~ l d ll n x n is the connection introduced in [I3|. Relation (J3) is a change of 
variables in the space of connections. Indeed, the original variables A M have 12 independent 
scalar functions. There are two independent scalar functions in n, four in C^, and, hence, 
there must be six independent scalar functions in W M . This is the case indeed because 12 
components of W ^ satisfy six independent conditions: Four in n ■ W M = and two in (||). 
Note that if is perpendicular to n, then covariant derivatives of with respect the 
connection a M + nC M is always perpendicular to n. The inverse transformation can be found 
by multiplying (^j) by n using first the dot and then cross products. The obtained relations 
allow one to express C M and as functions of and n. Substituting them into @, an 
equation for n as a functional of A^ is derived. It has been shown that for a given A M , the 
corresponding spin field can be computed as Q 

n = |tr (tQ^ a t 3 Q a ) , r = (n, r 2 , r 3 ) , (4) 

where 7$ are the Pauli matrices, tr(rjTj) = 25ij and Q A is a group element which depends 
on A^ so that the gauge transform of A^ with Q A satisfies the maximal Abelian gauge. 
Topological numbers of the defects have an integral representation via the spin field n 0. 

The gauge transformation Q A is in general singular (it might not even be single valued 
in spacetime). In other words, for a typical vacuum configuration of gauge fields, the max- 
imal Abelian gauge can only be achieved by a singular gauge transformation. According to 
numerical simulations J7|, the third component of the gauge fixed configuration A , which 
is associated with the unbroken £7(1) subgroup, carries Dirac magnetic monopoles. The 
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monopoles alone contribute about 90% to the energy of the flux tube (string) between static 
sources in the Yang-Mills theory. The monopole spacetime trajectories are determined by 
singularities of Qa- Hence, by taking Abelian connections of the monopoles and applying the 
inverse gauge transformation Vl\ to them, one can parameterize the relevant (or "monopole 
producing" ) degrees of freedom of a generic Yang-Mills connection. The result of this proce- 
dure is given by Eqs. (0) and (|3[). By construction, the spin field n carries all the information 
about the spacetime distribution of the defects, and, hence, its effective theory should describe 
the infrared physics of the original Yang-Mills theory. 

Using the path integral representation of the spin field effective action 0, one- loop calcu- 
lations have been done in [14j p|. They show that the action ([!]) is to be modified by adding 
terms containing higher order time derivatives (which might be a source of instability of knot 
solitons). The one- loop renormalization group flow for the effective action parameters also 
indicates that the effective action might be stable in the infrared region The results 



of [Tj] seems encouraging and deserve further studies, particularly, by some nonperturbative 
methods. 

Since the topological defects in the gauge-fixed theory are not solutions to the classical 
equations motion with a finite energy or action, it is rather hard to give them an accurate 
mathematical meaning in continuum quantum theory. The nonlinear change of variables 
(0) makes sense for classical connections whose values are well defined almost everywhere 
in spacetime. Typical quantum fields that form support of the path integral measure are 
distributions rather than classical smooth functions. The change of variables (0) is, in fact, 
ill-defined (because it involves products of distributions) unless some short-distance regulariza- 
tion is implemented. This can be achieved either by defining the path integral perturbatively 
with an ultraviolet cutoff, or by using the lattice (nonperturbatively defined) path integral 
for gauge theories. This latter approach is adopted in the paper. 

Here we develop the idea, first suggested in of using the inverse Monte Carlo method 
[16] to find out whether ([I]) is indeed a good approximation to the infrared Yang- Mills theory. 
Within the framework of lattice Yang- Mills theory, an explicit numerical algorithm is proposed 
to compute and study the effective action for the spin field and the renormalization group 
flow of its parameters, which comprises the main goal of the paper. 

2. In lattice gauge theories, the dynamical variables are group elements u X(l = u\ e SU(2) 
associated with each link, where x enumerates lattice sites, and /x indicates the direction of 
the link from the site x. Let {ui} be a Wilson ensemble of link variables distributed with the 
Boltzman probability Z^e~ Sw with Sw being the Wilson action of link variables and Zw 
the normalization factor (the partition function). The first step of the proposed numerical 
simulations is to generate an ensemble of the spin field from the Wilson ensemble. The spin 
field in the decomposition @ is defined by Eq. (Q). Thus, for every configuration u\ one 
has to find a configuration of gauge group elements fl x {ui) such that the gauge transformed 
configuration 

= Vt x ui^l\ + (5) 



2 see also jT5| where different approaches, not related to the observations in lattice gauge theories, have 
been explored. 
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satisfies the maximal Abelian gauge. Here x + fi denotes the lattice site next to x in the 
direction \x. The group elements Q x (ui) can be found by maximizing the function | 17| 



r 3 (uf)r 3 uf 



(6) 



for each configuration u\ of the Wilson ensemble. The collection of group elements Q x at all 
lattice sites is regarded as variables, while u\ are just parameters. For every configuration 
Ui the function \ u can have many local maxima. This is an evidence of the Gribov problem 
in lattice gauge theories (see for a review JE^] and references therein). For every element of 
the Wilson ensemble ui, one should take fl x (u) at which \ u attains its absolute maximum. 
Finding an absolute maximum of \ u is a difficult, if not impossible, task in the numerical 
gauge fixing. The state-of-the-art extrapolation toward the global maximum of \ u can be 
found in ||. The ensemble of the spin field is then computed as n x {u\) = |tr (TQl.T 3 fl x ) . 

It is also possible to find the lattice version of the change of variables (0) and therefore 
to obtain a system of cubic equations whose solution defines the spin field components as 
functions of link variables (the lattice analog of the equations for the spin field suggested 
in 0). Define an algebra element at each site n x = n\,T 3 Q x which satisfies the constraint 
trn^ = 2. Combining @ and (|6|) and introducing the Lagrange multiplier £ x to take into 
account the constraint on n x , the extreme value problem for (^) is equivalent to the extreme 
value problem for the function 

Xu(n) = Ex {E M tr [n x+f ,ujn x ui + (|trr^ - ljj . (7) 

Setting the variations of Xu{ n ) with respect to n x and £ x . to zero, the following equation for 
n, can be deduced 



ip x [n,u) 



~\~ ^x,i/ l ^x+^i'u x ^j 4- 2^77,3; . 



Now assume that fl x (u) is a local maximum of (§) and, hence, n x is a solution to then 
uy satisfies the maximal Abelian gauge and, by construction, n x = n x ■ r. The Lagrange 
multiplier £ x = £ x (u,n) is fixed by multiplying (||) by n x and taking the trace. After the 
substitution of £ x = ^ x (u,n) into fl3|), one gets the equation for the components of the spin 
field. Only two scalar equations in @ are independent. They determine two independent 
components of the spin field as functions of the link variables. 

As has been mentioned above, the group elements fl x and hence the spin field are not 
regular everywhere in space in the continuum theory. It is not difficult to give an example 
of Q x and the corresponding spin field n x such that the connection coincides with the 
Wu-Yang monopole (e.g., take n = x/r, r = \x\) which, in the maximal Abelian gauge, 
contains a Dirac monopole at the origin. The spin field is ill-defined at the position of the 
monopole. In the lattice gauge theory, the topological defects occur on the dual lattice sites 
|T7j| . For every configuration ui, the configuration Q x (ui) and, hence, n x are well defined and 



contain all the information about locations of the defects (magnetic monopoles) on the dual 
lattice and their topological numbers (magnetic charges). Consider an elementary cube of the 
spatial lattice and spins on its vertices. Frankly speaking, with an isolated defect present at 
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the cube center, the spins are directed outward (or inward) the cube. In the confinement phase 
the monopole-antimonopole pairs (or monopole loops, when describing topological defects by 
their spacetime trajectories) are condensed (no isolated defects), therefore the above simple 
visual picture would not be valid. However, n x is still well defined at each lattice site and its 
dynamics can be studied. 

3. The configurations of gauge fields u\ are distributed with the Boltzman probability 
Z^e~ Sw . The spin field configurations obtained from the Wilson ensemble must also be 
distributed with some probability Z~ 1 e~ s where S is the unknown effective action of the spin 
field. The problem is therefore: Given an ensemble of n x , find the corresponding probability 
or the effective action S(n). 

Any correlator of the spin field can be computed by the Monte Carlo method since the 
ensemble is known 

(F(n)) n = Z; 1 f Dne- s F(n) = J- £ F(n) + O(M^) , (9) 
J M {n} 

where the sum is taken over the ensemble of the spin field, Z s = J Dne~ s is the partition 
function, Dn = Yl x dn x and the integration over a spin at each site implies the integration 
over a unit two-sphere. Parameterizing the spin vector by the spherical angles 

n x = (cos (j) x sin 9 X , sin <fi x sin 9 X , cos cj) x ) (10) 

we get dn x = d<fi x d6 x sin 9 X where <p x G [0, 2ir) and 9 X £ [0,7r]. The expectation value (§) is 
also realized as an expectation value with respect to the original Wilson ensemble. Note that 
Eq.(p|) defines the spin field as a function of link variables n x = n x (u). Hence, 



[F(n)) u = Zw 1 j Due- s ^F(n(u)) , (11) 



where Zw is the partition function for the Wilson action. In principle, this observation can 
be used to determine the effective action directly via the original Wilson ensemble. The idea 
is the same as in the continuum case 0. Define the function A(u, n) by the condition 

[DnA(u,n)H6(<p x (n,u)) = l, (12) 

J X 

where (p x = ^tr (r<f x ) (see (§)), which leads to A(u,n) = det(d(p x /driy). The parameteri- 
zation (HID must be used to compute the derivatives and also the identity ip x ■ n x = is to 
be taken into account when computing the determinant. In spacetime the matrix dip%/dn b y 
appears to be sparse because nonzero elements can only occur for y = x, x ± /i. Substituting 
the identity (0) into the integrand in fllT]), changing the order of integration and comparing 
it with ([|), it is not difficult to deduce that 

S(n) = -\n(A(u,n)l[5(cp x (n,u))) , (13) 

where n x is now held fixed in the average over the Wilson ensemble. 
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Although QI5D defines the spin field effective action as an expectation value of some function 
of the link variables, we are interested only in its behavior in the infrared region. In the 
continuum case this amounts to the so called gradient expansion of the (nonlocal) effective 
action. In the numerical approach, the action is sought in the form 



S ~ ^ XiSi 



(14) 



where Si are some specified functions of the spin field, whereas the coupling constants Aj are 
to be determined. Since we are interested to compare the effective action with ([I]), we first 
find all possible local and independent terms that are of the forth order in derivatives and 
might contribute to the gradient expansion of S. In addition to the two terms in ([!]) it is 
necessary to include d 2 n ■ d 2 n and [(d^n) 2 } 2 . The four terms are all independent Lorentz 
and isotopic invariant terms containing up to four derivative operators. In fact, there is one 
more invariant term which can be built by contracting the Lorentz tensor d^n x d u n with its 
dual (like a ^-term in the Yang-Mills action) . Higher order terms can be classified accordingly 
by contracting invariant irreducible tensors of the isotopic and Lorentz groups with isotopic 
tensor products of the spin field and its derivatives. Then the renormalization group flow of 
the constants Aj must be studied as high momentum components of the spin field are removed 
(integrated out). The renormalization group flow would show whether or not the effective 
action is stable (or, in other words, is a good approximation to fljjD) in the infrared limit, 
and thereby validate or invalidate the conjecture. 
We set 

Si = — £^ (15) 

where rrii is the number of spins involved into a local interaction Si X . Let s Xjfl = n x+IJ ,—j xl j,n x , 
where j Xfl = n x+fJi ■ n x is defined by the condition s Xj/J/ 
with the continuum theory d^n ■ n 
are written as 



■ n x = to make the correspondence 
0. Let a be a lattice spacing. The local spin interactions 



S 



lx 



'T[s 2 + s 2 



S2X 
S3X 



S x,fi 



+ 



1 2 



x—/j,,ij, 



^x,fi X S x v j -\- \S X _ V ^ X .S' ,. ; 



) 2 + («. 



s 



Ax 



1 2 



£( 



y X,fJ, 



^X—flyfJ, ) 



+ [x — * x' = X + fi] + [x — > x' = X + jl] . 



(16) 

(17) 
(18) 

(19) 



In the continuum limit, the action ( |l~5f) goes into (|l|) with the additional terms described above. 
The #-term has the same form as S^ x where instead of the sum of squares, the sum of the dot 
products of each vector and its Lorentz dual has to be taken. Local spin interactions giving 
rise to terms with higher powers of <9 M in the continuum limit can be constructed similarly by 
using the correspondence rule: <9 M n — > 8 Xifl , 2d il d u n — > s x+lx>v - s XjI/ + s x+v ^ - 8 Xtfl , etc. 



6 



Note that two terms in S\ x give the same contribution as the sum over x is taken in 
15| ) and therefore mi = 2. Similarly, each of three terms in Si X (i = 2,3,4) gives the same 
contribution to (|15|) and rrii = 3. The reason the equivalent term are given in Si X is that Six 
is the part of the action Si that contains all terms involving the spin n x at a fixed site x. 
This representation will be useful in what follows. 

4. Here we formulate the inverse Monte Carlo algorithm for computing Aj. The inverse 
Monte Carlo method is well known in studies of the real space renormalization group of spin 
systems pi. It has also been applied to compute an effective action for monopole currents in 
the maximal Abelian projection fl8| . The use of the spin field order parameter n x rather than 



the monopole current is more appealing because of several reasons (relations to the quantum 
soliton theory, similarities between strings and knot solitons) pointed out after ([[]). 

Let S x denote all terms in 5" that contain the spin n x at a fixed site x, S x = J2i Aj«Sj X . For 
every Sj(n) we construct a new function 



Si(n, A) = — E Z* 1 [ dn x e- s *S ix = £ S i: 



(20) 



where Z x = Z x {n) = J dn x e~ Sx . The bar in Six denotes an expectation value carried out with 
respect to the effective action (■■■)„ but calculated for only one spin, n x . The environment 
(i.e. neighboring spins) is held fixed. So, S ix depends only on the spins at the sites neighboring 
with x, i.e., on n x ± fl and n x ±2^. Taking the expectation value of Si, we find the identity 



(S i (n)) u = (S l (n,\)) u . 



(21) 



Using the Monte Carlo method (g), the l.h.s. of (J2TD can be computed, while the r.h.s. cannot. 
The integral over n x in (|20|) cannot be computed for given configurations of neighboring spins 
because the true values of Aj are not known. Had the coupling constants been known, the 
ordinary integral in ( p0|) could have been computed, for instance, numerically for any given 
n x±)1 and n x±2fl . 

Suppose some trial values Aj of the coupling constants are taken to compute the r.h.s. of 
(|20|). The equality 

(Si(n,\)) v = (Si(n,\)) v (22) 

holds if and only if Aj = A^. This is used to set up an iterative algorithm to find the true 
coupling constants. Equation (22") is regarded as a system of nonlinear equations where the 
l.h.s. is known (cf. (0)). It can be solved numerically by Newton's method or some of its 



alterations. For Aj ~ 



A, we have 



(Si) u - (Si(n,X)) u 



(A, - Aj) 



(23) 



A=A 



Using the definition (|20|), it is not difficult to show that 

d ^ 



q Q Q CO 



(24) 
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The function SjSi is defined by ( E0"D where Si X is replaced by Sj X Si X . The true values of the 
coupling constants are computed by the iterating procedure 

fc(A (n) ) = E^( A(B) )( A i n+1) - A i n) ). ( 25 ) 

3 

k(X^) = (^(n)) u -(^(n,A^)) u , (26) 

-^j'C^ ^) = (^iSj — SiSj) u \x = \(n) , (27) 

where A^ = Aj and — > Aj as n — > oo. 

The convergence depends on the choice of the trial constants Aj. If iterations take many 
cycles, statistical errors are likely to introduce instabilities in the solution. A similar problem 



was encountered in |19 |. The solution there was to find the iteration limit through a linear 



mapping of the space of A*. Another method to compute the coupling constants is to use the 



Schwinger-Dyson equations ||20|| . In principle, the coupling constants can be compared with 
their "exact" values defined through (|13|). The expectation value in the r.h.s. of ([U|) can be 
expanded into a series over the spin field around some specific spin field configuration. The 
expansion coefficients can be computed my the Monte Carlo method in the Wilson ensemble. 
For instance, the mass scale Ai can be obtained by taking the second derivative of fll3| ) with 
respect to the spin field at the particular configuration n a x = 5 3a (as was suggested in the 
continuum case 0]). This procedure involves, however, computations of the determinant, 



which is very costly. Eq. (11) can be used to measure the goodness of the approximation 

dHH0). 

5. Having found the coupling constants, the renormalization group flow for them has to 
be investigated to prove the stability of the effective action in the limit of large wave lengths. 
With this purpose, we use the representation (|TTj) to take into account the constraint n 2 x = 1. 
Let the matrix f xk be a discrete Fourier transform associated with the lattice, J2x fk'xfxk — &k'k 
and J2k fx'kfkx = Sx'x- The sum over k implies the sum over all momentum vectors allowed 
by the lattice. Given the ensemble of 9 X and <p x , the Fourier components 9k and <pk can be 
computed. 

Next the spin field ensemble can be generated for all momenta bounded from above by 
some scale Ai 

^(AO = UkB k ■ (28) 
keKi 

Similarly for (p x (Ai). The sum in ( |2"8D is extended over those vectors k whose norm does 
not exceed the scale Ai. This subset in the momentum space is denoted K\. The ensemble 
n x (Ai) can be used as the input for the inverse Monte-Carlo procedure described in the 
previous section to compute new coupling constants A,(Ai). 

Repeating this procedure for successively smaller scales A^ +1 < A^ we can obtain the 
sequence of the coupling constants Aj(A fc ), k = 0,1,2,..., where k = corresponds to the 
coupling constants computed with the original spin field ensembles. By truncating the sum 
over momenta in (ESI) we generate a spin field ensemble in the infrared region (large wave 
lengths). Hence the sequence Aj(A&) describes the behavior of the coupling constants as 
functions of the scale A that restricts allowed momenta from above in the effective theory, 
i.e., as A decreases, the infrared limit is approached. 
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The explicit removal of Fourier modes can be strongly affected by the breaking of rotational 



symmetry on coarse lattices or for larger A's. So, the block spin decomposition JIB], |2]J might 
be a more attractive procedure to study the renormalization group flow. The idea is to average 
spins over elementary cells (blocks) of the original lattice. For instance, the angular variables 
9 X and <p x are specified at elementary cubic cell vertices. Consider a new lattice with mesh 
2a which is constructed as follows. Let the point y be the center of the elementary cell. The 
neighboring sites are then y ± 2\i. So, each site y of the new lattice is inside an elementary 
cube C y of the original lattice, and the cubes C y and C y i do not have common vertices if 
y 7^ y', while C y and C y > coincide if y = y' . Define 

6 y = 2-°- 1 J2 0, (29) 
xeCy 

and similarly for <f) y , where D is the lattice dimension. That is, 9 y is an average value of 9 over 
all vertices of one elementary cube. The spin field n y is defined by ([lOj) where x — > y. The 



averaging (^) is also equivalent to removing short wave length components of the spin field. 
Doing this procedure for a successively larger lattice spacing (2a, 4a, etc.) and computing 
the coupling constants on each step, we can again generate the renormalization group flow 
Ai(A fc ) (where A k ~ 2~ fc /a). 

The behavior of A, (A) allows one to verify whether the effective action (p]) (possibly with 



extra terms) is stable in the infrared region as was observed in [T%] in the one- loop approx- 
imation. For instance, it is critical to observe the right signs of Ai j3 (cf. ([!])) because Ai 
sets the mass scale for knot solitons, while A3 should reproduce the running gauge coupling 
constant g in the continuum limit. Relatively large (and growing) values of A2,4 would mean 
instability of knot solitons. It would also indicate that the approximation (|I]) is not justified 
and the higher order spin interactions are relevant for the spin field dynamics. In this case 
more higher order terms have to be included into (|T^) and the renormalization group flow of 
the corresponding coupling constants is to be computed. 

It should also be noted that the lattice size becomes important in studies of Aj(A). It is 
clear that A max = A ~ 1/a and A. min ~ 1/L where a is a lattice mesh and L = Na with N 
being the number of lattice sites in one direction. Values of A cannot be taken too close to 
A m j„ because the characteristic scaling behavior of the dimensional Aj in the continuum limit 
must still be observable. 



Let us summarize the essential steps of the algorithm. First, a Wilson ensemble of the link 
variables is generated. Then the spin field ensemble is computed. For every element of the 
ensemble the functions Si are computed for some trial coupling constants. The integral over a 
spin field at a fixed site x involved in the definition of Si has to be done numerically. Then the 
iterating algorithm is applied to compute the coupling constants. Finally, the procedure is 
repeated for several ensembles of the spin field which are obtained from the original ensemble 
by truncating (integrating out) short wave length components. This gives the renormalization 
group flow of the coupling constants which can tell us about the stability of the effective 
action in the infrared limit. Using an appropriate scaling of Ai in the continuum limit, one 
can compute the mass scale of knot solitons supported by the effective action (|T|), while the 
coefficients at the additional terms would determine their stability. 
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As a final remark, it should be noted that the choice of the spin field ensemble suggested 
above is associated with the maximal Abelian gauge. It might be interesting to investigate 
the effective action if different Abelian gauges are used to obtain the spin field ensemble. 
The purpose of such a study is to seek a numerical evidence that the effective dynamics 
of Yang-Mills fields in the infrared region is indeed governed by configurations that exhibit 
topological defects when taken in any Abelian gauge |23| . A generalization of (|S|) for the spin 
field ensemble can be defined as follows. Let Q x (u) be a gauge transformation such that -up 
satisfies some Abelian gauge. The latter implies in particular that there exists a field defined 
at each x, $ x = <& x {u n ) which is diagonal [r 3 ,<iy = and also <& x {u n ) = Vt x § x {n,u)VL x . 
Therefore, a general equation for the spin field reads 

[$ x (n,u),n x ] = , (30) 

where Q x is any function of n x and ui which transforms according to the adjoint representation 
of the gauge group. An example of the field Q x (n, u) is given by the first two terms in Eq. (^), 
and equation (§) itself is clearly equivalent to (|30D . Link variables in two different Abelian 
gauges are related to each other by a singular gauge transformation. Equation ( |30"D is covariant 
under such transformations and, hence, represent the most general equation for the spin field 
components. 

The simplest class of Abelian gauges is obtained by taking <& x to be independent of the spin 
field. In this case, the spin field is identified as a solution of (|3U[): n x = § x (u)/\$ x (u)\ where 
the vertical bars mean the norm with respect the trace scalar product in the Lie algebra. This 
is the case of the Polyakov Abelian gauge |24| and the so called Laplacian Abelian gauge |25 



The spin field ensemble in the case of the Polyakov Abelian gauge is probably the simplest one 
to generate. In fact it is possible to average over all such choices of $ in the Yang-Mills path 
integral (dynamical Abelian gauge |26[) by using a supersymmetric extension of the theory. 



The SU(N) generalization of (pOl) is also straightforward. Take N — 1 orthonormal and 



commutative spin fields, (n x ,n x ) = 5 ab and [n",n^] = where a,b = 1,2, N — 1. They 
should satisfy the equation [3> x (n, u), n x ] = where $ x is some operator that transforms 
under the adjoint representation. In the simplest case when <& x depends only on the link 
variables, the solution for the spin fields is obtained by orthogonalizing the commutative 
fields ~ dhl-ia+^K 2 • --K a+1 where e i is an orthogonal basis in the Lie algebra su(N), 
$ x = e l Q l x and d!fj 2 ... ia+1 are the N — 1 symmetric irreducible invariant tensors of su(N). Since 
for any Lie algebra element $ x it always possible to construct N — 1 linearly independent 
elements (using the tensors that commute with <& x and amongst each other, the spin 
field can also be specified by [Q x (n,u),n x ] = for some choice of & x (n, u). A generalization 
of Eq. (^|) to the SU(N) case is simple. One has to replace n x by n x in the first two terms 
(this defines $"(n, u)) and £ x n x by LbCf^ in the third one where £f = (^> x ,n b x ). 
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Note added. After submission of the paper for publication, I have learned that the group 
of A. Wipf have done numerical simulations on a 16 4 lattice The mass gap in the spin 



field spectrum and possible global SO (3) symmetry breaking (a preferable direction of the 
spin field) have been reported so far. The existence of the mass gap around the lowest glueball 
mass should certainly be expected because of the Abelian (and monopole) dominance. The 
interesting question of the renormalization group flow of the coupling constants is still to be 
studied. 
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